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ABSTRACT 


In modern turbofan engines, variable geometry has been 
incorporated to improve some off-design performance. Most 
control designs ignore this variable geometry and use fuel 
metering as the primary control input. 

This thesis investigates the use of variable geometry to 
control the engine and, thereby, reduce fuel consumption due 
to transients. Additionally, steady-state trim conditions 
are altered to reduce the static fuel consumption. The non- 
linear transient simulation program is used to analyze the 
steady-state operating condition and develop small perturba- 
tion control limitations. Linear models, both large and 
reduced order, are used in analyzing the effect of controllers 
on system response. A computer program was generated to re~- 
duce a large order linear model to a usable size for control 
system development. 

This analysis shows the reduced-order regime dependent 
controllers to be viable and to favorably enhance the quest 


for reducing specific fuel consumption in existing engines. 
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I. INTRODUCTION 


High performance, military aircraft are operationally 
tasked to conduct missions throughout their designed per- 
formance envelopes. To accomplish these missions, the air- 
craft must operate efficiently within its flight regime. 

The typical jet engine is designed for operation at a single 
specified altitude and Mach number combination, e.g., 30,000 
feet at M= 0.9. This design point will give inherently 
good cruise performance and is particularly well-suited to 
commercial transport operations. This engine design method 
forces the user into accepting off-design performance in 
typical military applications. Engine performance needs to 
be optimized to efficiently conduct all assigned missions, 
e.g., air combat maneuvering, ground attack, supersonic inter- 
ception and carrier operations. Adaptive digital control 
methods can be incorporated within the engine to accomplish 
this mission dependent optimization. 

Modern control methods allow the designer to develop a 
control system to regulate every parameter of the parent 
system. He has been aided by computer simulations of non- 
linear systems and linear approximations of those systems. 
The desired control system can be developed parallel to or 
ahead of the controlled system. 

The state-of-the-art turbofan engine utilizes variable 


geometry, in such forms as variable incidence inlet guide 


ll 


vanes, adjustable exhaust nozzle area and airflow bleed-off. 
Hydromechanical controls have been used in the past to im- 
prove some elements of engine performance through positioning 
these variable engine components. Electronic controls, such 
as were introduced in the Pratt & Whitney F100-PW-100 turbo- 
fan engine, allow improved scheduling of the variable devices, 
again improving the engine performance. Studies have been 
done on the F100 engine to optimize its performance at a 
point other than its design point [{1l]. 

Originally, the onboard microprocessing capability of an 
aircraft was intended for weapon system improvement and 
integration. Further computing equipment was added as part 
of the engine control to supervise the hydromechanical con- 
trol system and detect faults. The F100 engine introduced an 
electronic control which allowed use of digital control tech- 
niques incorporating the hydromechanical devices as implemen- 
tors rather than controllers. Digital device improvements 
continue to expand the available computing power available, 
making flight control and engine control integration a dis- 
tinct possibility. As more software space becomes available 
to the engine controller, regime adaptive control will be the 
accepted standard in military applications. 

This thesis looks at using modern, multivariate control 
logic in designing flight-regime dependent control for a 
typical military turbofan engine, the Pratt & Whitney F100. 
Specifically, the attempt is made to improve thrust specific 


fuel consumption, i.e., fuel economy, in the high altitude 
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cruise environment. Such improvement in fuel economy is 
especially important in the conservation-conscious decades 


ahead. 
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II. TURBOFAN ENGINE PERFORMANCE CONSIDERATIONS 


A. FLIGHT REGIME FACTORS 

The ability of a military aircraft to perform an assigned 
mission is determined by its flight regime. Figure II-l is 
a typical flight regime for a modern, high-performance 
military aircraft. 

The confines of the flight regime are dictated by numer- 
ous aerodynamic and propulsive factors. A well designed 
aircraft will incorporate an engine having operating boundar- 
\ ies in excess of those of the airframe. This design technique 


allows for a factor of safety in actual operation as well as 


compensate for installation and equipment aging losses. 

At the low speed, low altitude region of the envelope, 
engine performance is judged in the ability to produce 
aufficient thrust for take-off and landing and to provide 
rapid engine acceleration in the missed approach or balked 
landing. The extreme upper region of the flight envelope 
requires the engine to produce maximum thrust while not ex- 
ceeding turbine temperature limitations. 

Figure II-2 is a mission delineated flight regime for the 
supersonic tactical fighter aircraft. The desired mission 
performance from this chart must be interrelated with the 
aerodynamic and propulsion limitations of Figure II-l. An 
adaptive control system can be designed to schedule optimum 


performance of the airframe and engine throughout the flight 
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regime. Northrup Corporation is involved with a U.S. Air 
Force project to implement an adaptive flight control system 
in an aircraft such as the F-16. The same type of approach 


can be applied to engine control and integrated into an 


aircraft control system to optimize the total package. 


B. TURBOFAN ENGINE SRFORMANCE MEASUREMENTS 

The primary factor with which to judge engine performance 
is the engine's ability to convert fuel to thrust efficiently. 
This measurement, thrust specific fuel consumption, SFC, is 
a function of the airflow through the engine, Mae the bypass 
ratio, b, and the fan and core exhaust velocities, wo and | 


Ue": also termed cold and hot stream velocities, 
M, 
SFC = pl (i+f) 0 "+ 3U_'' = (14+8)V]) (II-1) 
£ e e 


where Me is the fuel flow, £ is the fuel-to-air ratio and V 
is the flight velocity. Figure II-3 shows the F100-Pw-100 
engine in cross section with reference labels. 

The engine designer fixes the airflow and bypass ratio of 
the engine. The exhaust velocities are functions of the ex- 
haust temperatures and pressures, which are controllable 


through fuel metering and variable geometry. 
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Figure II-3. Engine Cutaway 


where y and ¥, are the respective ratios of specific heats, 
To6 and Tog the respective exhaust total temperatures, is 
the exhaust nozzle efficiency, and Pa/Po6 and P./Pog the 
respective ratios of free stream ambient pressure to dis- 
charge total pressure. 

The standard bypass fan will produce a fixed pressure 
ratio and exhaust temperature, producing an exhaust velocity 
in the cold stream in excess «* the flight velocity. This 
increases the engine's thrust without altering the hot stream 
through the engine core (ignoring energy lost in driving the 
fan). The advanced, low-bypass fans that reintroduce the 
cold flow into the augmentor, as in the F100, increase the 
augmentor efficiency and during non-augmented operation, 
allow for increasing ue" through the use of the variable ex- 
haust nozzle. 

Studies show that Thrust Specific Fuel Consumption, SFC, 
can be decreased through increased turbine inlet temperature 
and compressor pressure ratio. Both of these variables are 
designed into the engine at the engine design point. Engine 
control hardware can be designed to maintain the turbine inlet 
temperature, FTIT, and pressure ratio, Ro at the optimum or 


best achievable level during off-design operations. 


C. ENGINE CONTROL DEVICES 
The mechanical devices available for engine control are 
the hydromechanical fuel control, the variable exhaust nozzle 


and variable incidence inlet guide vanes or stators. The 
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first device has the effect of increasing the fuel-to-air 
ratio and temperatures in the turbine and exhaust sections of 
the engine. The combined effect of increasing f and ees is 
increased thrust, but does not necessarily improve SFC due to 
the corresponding increase in the fuel flow. The variable 
area exhaust nozzle has been used for the last two decades in 
most high performance military aircraft. This device in- 
creases the hot sectior ei. cst velocity, Par thus increasing 
thrust and maintaining ix: »‘@al pressure in the exhaust 
section of the engin. 

The latest mechanical devices are the variable incidence 
blading concepts. fhe tirst use of these devices was to 
achieve near design mass flow rate through the compressor 
during the engine start cycle. These variable stators would 
be rotated to a starting configuration until the engine reached 
a given rotational speed and then be driven to a fixed normal 
operating position. 

The full advantage of variable geometry blading comes in 
the optimizing of off-design performance. The single stage 


pressure ratio is found to be 


2 


Po3 ue Cc, 
; — = l+ny ape! - g (tan b, + tan a,)] (II=4) 
| Pou ol r 


where » is the stage efficiency, U. is the rotational velocity 
of the blading, R is the Universal Gas Constant and To. the 


total temperature at the front face of the stage. At the 


Scr 


design condition, the factor (tan by + tan ay) is minimized, 
thus providing the best pressure ratio. The angle bo is the 
relative flow direction upon leaving the rotor portion of the 
machine; ay is the angle at which the flow impinges the rotor 
and is a function of the stator angle (see Fig. II-4). Con- 
trolling the stator angle increases the stage pressure ratio 
in the early stages and increases the compressor (or fan) 
performance over a wide range of operations. 

In the F100-PW-100 engine, four mechanical devices are 
used, fuel control, variable exhaust nozzle, inlet guide vanes 
with moveable trailing edge and moveable stators in the first 
three compressor stages. The exhaust nozzle maintains nozzle 
area, expansion ratio and boattail drag, simultaneously, near 


optimum. The inlet guide vanes are used ahead of the fan to 


improve inlet distortion tolerance, improve fan efficiency 
and enhance engine acceleration performance. The variable 


stators improve starting and high Mach number characteristics. 


D. ENGINE MODELING 

It is not economically practical to build an engine and 
: then conduct experimentation to determine its performance. 
Modern technology, in both digital and hybrid computers, 
allows the design engineer to build a simulation model of the 
engine and predict the system response well in advance of 
hardware manufacture. Several generic computer simulations 
are in use, these allow non-specific engine performance analy- 
sis. Additionally, a computer model is generated for each 


‘ new development engine by the engine manufacturer. 
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Figure II-4. 
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Szuch, [2], provides a detailed description of current 
engine modeling techniques. The Pratt & Whitney F100-PW-100 
engine has been modeled with both a real-time hybrid computer 
[3], and a digital transient simulation [4]. Both programs 
have distinct advantages and applications. Although the digi- 
tal program is not real-time, it provides a large volume of 
data and is easily accessed and maintained. A copy of the 
F100-PW~100 Transient Engine Simulation Deck (CCD 1103.3.0) 
was obtained from the YF100 Special Project Office of the 
Aeronautical Systems Division of the U.S. Air Force at Wright- 
Patterson Air Force Base, Ohio. This program was modified 
to be compatible with the Naval Postgraduate School IBM 
370/3033 computer system. 

Designing a controller based on the non-linear digital 
simulation alone is extremely difficult, if not impossible. 
Miller and Hackney [5], developed high order linear approxi- 
Mations, at steady state operating points, from the non- 
linear simulation program. The linear system was modeled in 


the classical state variable form 
‘ X = FX + GU (II=-5) 
Y = HX + DU (II-6) 


The 16 system states, listed in Table II-l, were chosen to 


coincide with the states measured in the non-linear simulation. 


These states include 2 rotational speeds, 3 internal total 


—_ HY 


TABLE II-1 


Engine States 


xl -- Fan speed, (Nl), in RPM 

x2 <-~ Compressor speed, (N2), in RPM 

x3 -~ Compressor discharge pressure, (Pt3), in psia 
x4 -~ Interturbine volume pressure, (Pt4.5), in psia 
x5 -- Augmentor pressure, (Pt7m), in psia 


x6 -- Fan inside diameter discharge temperature, (Tt2.5h), 
in Rankine 


x7 <-- Duct temperature, (Tt2.5c), in R 

x8 - ##Canpressor discharge temperature, (Tt3), in R 

x9 -- Burner exit fast response temperature, (Tt4hi), in R 
x10 -- Burner exit slow response temperature, (Tt4lo), in R 
xll -- Burnerexit total temperature, (Tt4), in R 


x12 -- Fan turbine inlet fast response temperature, (Tt4.5hi), 
in R 


x13 -- Fan turbine inlet slow response temperature, (Tt4.5lo), 
in R 


x14 -- Fan turbine exit temperature, (Tt5), in R 


i x15 -- Duct exit temperature, (Tt6c), in R 


x16 -- Duct exit temperature, (Tt7m), in R 


pressures and 11 total temperatures. The inputs to the sys- 
tem are the fuel flow, the 3 variable geometry elements pre- 
viously discussed, and the customer bleed air percentage. 

The outputs include the net thrust, the engine total airflow, 
turbine inlet temperature and the stall margins. Addition- 
ally, two fan exit pressure ratios are determined for some 
of the operating points. The inputs and outputs are listed 
in Table II-2. 

The technique used by Miller and Hackney perturbed each 
state, xi, slightly while holding the other states and inputs 
constant. This allows calculation of the deviations caused 
by that state in the other states and outputs. These devia- 
tions then form the F and H matrices of equations (II-5) 
and (II-6). Perturbing each input, ui, in the steady state 
configuration determine the G and D matrices. The sampling 
time was set at 7 milliseconds. 

The linear models obtained are very good approximations 
of the non-linear system but are often too complex and/or do 
not contain the most convenient parameterization to be used 
in a practical design. A reduction of model order and, per- 
Raps, augmentation of the linear model is generally required 


prior to using the model for control design. 


E. CONTROL DEVICE EFFECTIVENESS 
Before designing an engine control, one must determine 
the effect each input has on the system and the responsive- 


ness of each input. The five control inputs listed in Table 
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TABLE II-2 


Engine Inputs 


Main burner fuel flow, lb/hr 


Nozzle jet area, ft**2 


Inlet guide vane position, deg 


High compressor variable vane position, deg 


Customer bleed flow, % 


Engine Outputs 


Net thrust, lb 


Total engine airflow, lb/sec 


Turbine inlet temperature, R 


Fan stall margin 


Compressor stall margin 


Fan exit delta p ratio (test data) 


Fan exit delta p ratio (theory) 


II-2 are used in controlling the F100 engine. The customer 
bleed air percentage, u5, is essentially constant in steady 
state operations and can therefore be ignored in the control 
development. A study could be made to determine the merit 
in extracting additional bleed air from the engine as a pri- 
Mary control, and dumping that portion not required by the 
customer systems. 

The intent of this control development is to minimize the 
specific fuel consumption of the engine. At the same time, 
only small variations in turbine inlet temperature, y3, and 
total engine airflow, y2, are desirable and no change in 
thrust is allowable. One method might employ using exhaust 
nozzle area, u2, and fan inlet guide vanes, u3, inputs to 
minimize SFC at the trim condition and then compensate for 
thrust loss by using fuel to balance the equation. 

Figure II-5 plots the change in specific fuel consumption 
due to incremental changes shown in Table II-3 in the inputs 
ul, u2, and u3; input u4, the compressor variable guide vanes, 
exhibit no influence on SFC. Proportional derivatives, i.e., 
slopes were determined by making a linear approximation to 
the curves. It is clear in the figure that nozzle area and 
inlet guide vane inputs (u2 and u3) have the most substantial 
effect on SFC. 

Figure II-6 plots the net engine thrust, yl, versus the 
control input changes. Again u2 has the greatest effect and 
ul and u3 contribute about equally to thrust variations. It 


is also noted that a linear combination of ul and u3 can be 


aS ek 


syndu]I To12uU0D SA O4S UT sbueYyd °G-II ainbtd 


gnduzr [or,uUO0D 


0 


ee ae a er 4 ---—--- —~- —--. ---$ - ‘i 


p/en 


+ ene ee eM ee ee  --- 


Seiiiensesienttn eacmanemeeaneeatiees -- ie 


———— 


O4SV 


T0° 


28 


Syndur JTOI}UCD SA ysNaAYyL UT ehueyd °9-II sanbtg 


qndujy fo2z}u0D 


[tA] (aT) 
ASnNIAQLV 


29 


00zc 


syndu[ [To1}u0D SA 
ainjzeraduay, WeTUT sutTqing, ut sHueyD 


yndul jyor3uU0D 


- eeheiedieeheeattieed sitertnediieemeneatiaaens | eelis eter REET t ae tee 


"L-II eanbtga 


i} 


(¢ 4) 

(De) 
91ny 
-eraduaL 
yaptuyl 
auTqINLV 


30 


TABLE II-3 
Input Perturbations 


The following perturbations are used throughout this 
investigation. They represent the largest expected value, 


or physical limitation placed on 


each control input. 


Fuel Flow, ul + 200 lb/hr 

Nozzle Area, u2 + 0.2 sq. ft. 

Inlet Guide Vanes, u3 + 4 degrees 

Compressor Vanes, u4 + 4 degrees 

Bleed Air, u5 0.0 percent 
used to offset the thrust change caused by a change in u2. 
This is important in meeting the zero net thrust change de- 
sign requirement. Figure II-7 shows the effect the inputs 
have on turbine inlet temperature, y3. Again, u2 is the 


dominant control input. These figures clearly demonstrate 
the need for minimizing the exhaust pressure mismatch with 
ambient conditions. 

Table II-4 is a summary of the proportional derivatives 
extracted from the preceding figures. Using these derivatives 


the following relationships are found. 


G(SFC) = 2.27E-05 d(ul) + .0655d(u2) 
- .00092d(u3) (II-7) 
A(yl) = .9325 d(ul) - 1351.7d(u2) 
+ 18.6275d(u3) (II-8) 
31 


TABLE II-4 
Summary of Derivatives 


Input SFC Thrust (yl) FTIT (y3) 
Fuel Flow, ul 2.27E=5 -9325 12126 
Nozzle Area, u2 -0655 -1351.7 283.085 
Inlet Guide Vanes, u3 -.00092 18.6275 -2.1641 


U4 was found to have no effect on the parameters of 
interest. 


d(y3) = .12126d(ul) + 283.085d(u2) - 2.1641d(u3) (II-9) 


If a decrease in u2 of -.2 square feet and an increase of +4 


degrees in u3 are applied, a change of -.017 occurs in SFC. 


The fuel flow, ul, is used to balance the thrust equation. 

To accomplish this the fuel is reduced by 370 pounds per hour, 
further decreasing specific fuel consumption. The total effect 
of these three inputs at the trim condition is a 3.5% reduc- 
tion in SFC, no thrust change and a 110 degree decrease in 

the turbine inlet temperature. 

By requiring that thrust level not change, the SFC per- 
formance is compared under identical operating conditions. 
Each input can be studied separately to determine its effect 
on the outputs and on SFC. To remain in the small perturba- 
tion regime inputs should not be varied more than 10% of their 
steady-state values. Combinations of inputs are then chosen 
to effect the most desirable improvement in SFC and minimize 


the effect on engine performance. 
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This technique has shown improvement in the specific fuel 
consumption at the operating point chosen in the studies of 
Reference 1. The use of the existing variable geometry to 
improve the fuel efficiency at the high altitude cruise 
operating point chosen for this investigation will be dis- 


cussed in a later section. 
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CONTROL DEVELOPMENT 


A. LINEAR OPTIMAL CONTROL THEORY 


Given the linear system, 
X = FX + GU (III-1) 
consider a linear control law defined by 
U = -KX + Uy (III=-2) 
The original system can be modified to 


X = (F - Gk)X + GU, (III-3) 


The response of the modified system, as well as its eiyen- 
values, is determined by the effect of the gain matrix, K. 
For a low order system, the control gains could be found by 
trial-and-error methods and applied to the system to determine 
if the response was acceptable. This method is impractical 
for large order systems. 

Optimal control is achieved by minimizing a cost function, 


J, which is defined as 
J = $f [(X*R,X) + (U7R,U) }at (III-4) 


where Ry and Ro are symmetric weighing matrices; Ry is non- 


negative definite and Ro is positive definite. A third portion 


can be added to (III-4) to provide terminal state weighing, 


T 
Xe, PK (ty) (III-5) 


where Py is a non-negative definite symmetric matrix. 
It is shown by Kwakernaak and Sivan [8], that the optimal 
input, Ug (te), based on the cost function is given by 


Uy(t) = (R,)7*e™P(e) (111-6) 


thus the knowledge of p(t) solves the regulator problem. A 
2n linear system is formed with the optimal system behavior, 


X(t), and the adjoint variable, p(t). 


: | 
-1_T 
Xp (t) F -GR,G XQ (t) 
= (III~7) 
p(t) -R, -F? p(t) 
where 
p(t) = P(t)X, (t) (III-8) 
with 
-1 T 
K(t) = Ro (t)G™ (t) P(t) (III=9) 


P(t) is a symmetric non-negative definite matrix that 


satisfies the Riccati equation 


-P(t) = Ry (t) - P(t)G(t)R>*(t)G" (t) P(t) 


+ P(t) F(t) + Po (t)P(t) (IIr-10) 


with the terminal condition of 


P(t,) = Py (III-11) 
Under these conditions, as the terminal time, tl, approaches 
infinity, the control law asymptotically approaches a unique, 
Stable steady state condition. Thus, the matrix Riccati 


equation, (III-10), can be rewritten 


1 T. 


Oo = Rk -pcriic™p + pr + F Tp (III-12) 


2 


and 
K = R-“G'P (III-13) 


Many numerical methods are available to solve the matrix 
Riccati equation for the steady-state control gains. Some 
of these are covered by Kwakernaak and Sivan [6]. This 


thesis employed a computer solution, OPTSYS4, developed by 


Hall and Bryson [7], at Stanford University. A representative 


output is included in Appendix B. 


B. CONTROL SYSTEM MODELS 
At this point, an analysis of control models is made to 


establish the form of the controller to be applied. A 


- ee ee ee 


a 


tet OM 


comparison of small and large perturbation controls must be 


made. Take the linear system 


X = FX + GU (III-14) 


K 
u 


HX + DU (III~15) 


This system can be made closed-loop by incorporation of a 


control law 
U = = KX (III ~16) 


K can be a time-varying matrix or constant as required by 


the system. A more general control law would be 


Uo = UL + C(x =x) + Jcyacy ~Yp)dt (III-17) 


where A is the diagonal matrix whose diagonal terms are l 
or 0 determining which outputs are used in formulating the 
input. Cy. and cy are again matrices associated with the 


states and the outputs, respectively. Additionally define 


U = U, + Rit (III-18) 


“ 
) 
“s 
+ 
ye) 
“ 
ct 


(III-19) 


the 0 subscript indicates trim position/value; “a and R, are 


functions of u and x associated with time. For small pertur- 


bations the time varying terms, Ry and Ry! can be ignored and 


ie ee ee 


equation (III-17) can be rewritten 


su = CL 5X + J CyAsY dt (III-20) 


Dropping the delta terminology and taking the derivative of 


this and using equation (III-16), yields 


- KX = CX + CyAY (III~-21) 


combining terms 


(-K-c )X = Cay (III-22) 
x y 


Substituting, using the system description of equations 


(III-14) and (III-15) one obtains 


(~K-C.) (FX+GU) = CA CAX+DU) (III-23) 


and solving for Cy. 


cs “K-C A (H~DK) (F-GK) “> (III-24) 


and using all outputs (A = I) and equally weighing them 


(Cy = I) one has 


C. = -K-(H=DK) (F-GK)~* (III-25) 


However, if the outputs are ignored (A = 0), i.e., consider- 
ing the case where small changes are made in the output 


vector, then 
Cc = - K (III-26) 


as in the simplified control law of equation (III-16). 

This approach does not apply to the large perturbation 
case. In large perturbation analysis, the time varying 
terms cannot be ignored. The mean values of u and x are as 
in equations (ITI-18) and (III-<-19). R, and R, are the time 
functions that describe the path that the mean values follow 
over the period of the large perturbation. The control law 


becomes 
sU-Rt = C(x -R,t) + / CLASY dt (III-27) 


Again, if all outputs are used and equally weighed, then 


taking the derivative 
“RK “Rt = Cx =C RY “CR t een (III-28) 
and combining terms 


~C,, (F-GK) X +R, +RC Lt = K(F-GK)X + Ry 


+ Rt + (H-DK) X (III=29) 


If we set 


Rc = R (III-30) 


and 


RC = R (III-31) 


in effect, assuming a linear relationship between the states 
and inputs. Then 

-C, = -K -~ (H-=DK) (F-GK) 7? (III-32) 
as in equation (III-25) for the small perturbation case. 

An adaptive control system removes the need to mojel the 
large perturbation response. Linear models can be developed, 
as described in Section II.E, to define the large perturba- 
tion case as a series of linear relations. The onboard digi- 
tal computer then schedules the control gains based on the 
linear models as the flight conditions pass through the model 
transition points. The rapidity of response inherent to the 
digital computer removes the requirement for modeling the 
total dynamic system, rather, only the dominant response model 
for the selected operating point. The reduced-order linear 


model is far easier to handle in formulating the control sys- 


tem and is derived to match the dominant system characteristics. 


C. REDUCTION OF LINEAR MODELS 

DeHoff and Hall [8], present a method of reducing the 
order of the linear system based on the dominance of the 
states and the crosscoupling of the states as observed through 


eigenvector analysis. The required order of reduction is 


determined by the number of dominant states appearing in the 
desired control bandwidth. For the F100 study, DeHoff and 
Hall chose a control bandwidth of 0 to 10 rad/sec, corres- 
ponding to the primary control device (fuel flow actuator). 
Their operating point (static, sea-level, intermediate power) 
dictated that the control be designed to modulate thrust, 
providing maximum thrust for take-off performance. At another 
operating point a different parameter may dictate the control 
design and, thus, the primary control actuator and control 
bandwidth. 

The method of decomposition involves finding the eigen- 


vector transform matrix, T, such that 


XT = TA (III-33) 


J where \ is the diagonal eigenvalue matrix. Then, by defining 


an alternate state vector, Z, where 


X = T2 (III-34) 


the original system may be written in modal coordinates as 


(IIIT-35) 


If a control bandwidth is determined, 


the eigenvalues 


falling in that range are the only ones that need to be in- 


cluded in the initial reduced order model. 


The matrix may 


be reordered by constructing a matrix Q, which has a 1 in 


the position corresponding to the eigenvalue to be included 


(column) and the reordered position (row) of that eigenvalue. 


Thus, if eigenvalue 5 of the original A matrix was to be re~ 


ordered to the first position, a 1 would be placed in the 


(1,5) position of the Q matrix. Now, 


A' = QA 


is incorporated in the linear model 


ez = A'z + or tau 
letting 
Z' = QZ 
2' = itgutg: + or ty 
and 
x = TQ 42: 


the reordered matrix 


(III-36) 


(III=37) 


(III-38) 


(III ~39) 


(III~40) 


(III-41) 


SSE Ra ee ram. Sag Si, 


This X' vector is the reordered state vector, where the 

elements contained in the subvector Xr are those states 

associated with the eigenvalues in the control bandwidth 
and X, contains all other states. 


The linear system may be rewritten 


x! ~ RFR Lx! + RGU (IITI~42) 
Now define 
FP’ = RFR + (III-44) 
G' = RG (IIIT-45) 
H' = aR (III-46) 
and 
' 
es te es 
FP‘ = woresnon-- (III-47) 
t 
For « Fo2 
mk 
G' = ca eas atl ella (III-48) 
3 
H 
H = (Hy Ho] (III~49) 
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If the states in X) truly model the full order system, 


then Xo will be essentially zero and 


x) = Fiyky + Pi2%s + GU 


Solving for X, in equation III-51, one has 


2 


-1 


Ro =~ Foo Fou%1 * Fae 


This is incorporated into the first equation 


e 


7 7 zi R ai 
X= (Pi -FioFoo Foy) %1 + (G)-Fi2F22 G2)U 
and the output equation becomes 

Y = (H,-HoF.>'P..)x, + (D-H F,74G,)U 
p7HaFo2 P27) *y 2F22 S2 


The reduced order model can be written as 


Xx = Frx, + Gru 


BF 


where 


me eee 7 wig ot ies 


. 


(III-50) 


(III-51) 


(III-52) 


(III-53) 


(III-54) 


(III-55) 


(III=56) 


(Lit<=S7) 


= -= =1 - 
Gr G, - FyoFo3 65 (III-58) 
Hr = H_-HoF.7lF (III-59) 
17HoF 22 Fay 
Dr = D-H.F.746 (III-60) 
2F22 S2 


This technique was used to develop the program REDUCR2, 
documented in Appendix A. REDUCR2 was used to develop the 
reduced order models used in the control development portion 
of this thesis. Sample input and output of REDUCR2 are 


listed in Appendix A. 


D. REDUCED ORDER MODELS 

Three reduced order models were generated with the REDUCR2 
computer program (Appendix A). The criteria for selecting 
the desired states was that used by DeHoff and Hall, a con- 
trol bandwidth of 0-10 Hertz. The eigenvalues of the full 
16th order system were found using the OPTSYS4 computer pro- 
gram and then associated to the system states using modal 
analysis. Table III-l gives the results of this modal analy- 
sis. The five states associated with the control bandwidth 
are fan speed (xl), compressor speed (x2), augmentor pressure 
(x5), burner exit temperature (x10), and fan turbine inlet 
temperature (x15). 

A fifth order model was generated using these states. It 
was noted that only three of these states (the two rotational 


speeds and the pressure) were dynamically significant, i.e., 
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TABLE III-l 


Operating Point Eigenvalues 


State 

X4 P04.5 
X3 P03 

X11 T04 

X16 T07 

x9 TO4Hi 
X12 T04.5Hi 
X7 T02.5C 
X14 TOS 

X15 TO6C 

X6 T02.5H 
X8 T03 

x5 P07 

X1 Nl 

X13 T04.5L0 
X2 N2 

X10 NO4LO 


Physical Property 


-540.23 
-156.43 
-49.69 


- 54.12 
- 42.61 
- 22.59 


-21.58 


18.62 


14.49 
7.35 
3.60 
1.98 
1.31 

64 


Eigenvalue 


2 


+-313599 


4 92592 


control 
band- 
width 


7 


—s ~« 


TABLE III-2 


Reduced Model Comparisons 


Model #1 Model #3 Model #2 Full Order 
X21 173.51 175.89 175.89 173.41 
X2 125.10 126.86 126.86 125.02 
x5 . 326 «32° .325 . 326 
X10 4.85 4.86 
x11 53.51 53.48 
X13 3.98 3.98 
X16 23.04 23.03 
yl 94.48 174.8 97.70 156.99 
Y3 11.84 53.51 9.3 53.48 
_ 4? 


represent energy storage in the form of torque and pressure. 
A third order model was generated using these three states 

, and a comparison made. In this initial analysis these two 
models produced virtually identical response with respect 
to the key outputs thrust (yl) and turbine inlet temperature 
(y3), see Table III-2. The OPTSYS4 program was used to deter- 
Mine the transfer functions associated with these two outputs 
for each model. It was found that the fifth order model con- 
tained zeros near the poles at -1.960 and ~.651, the eigen- 
values associated with the two temperatures in the fifth 
order model. This caused these two poles to have no effect 
on the response of the model. Thus the third order system 
sufficiently models the control bandwidth. 

In both models, the turbine inlet temperature (y3) response 

did not closely model that of the full order system. There~ 
fore, a third model was formed using the three states of 


the third order model plus the conjugate pair of poles repre- 


senting the turbine inlet temperature (xll) and exit tempera~ 
ture (x16). This model gave excellent thrust and temperature 
response (yl and y3, Table III-2) due mainly to the inclusion 
of the turbine inlet temperature as a system state. Table 
III-3 lists the three models, their states and associated 
eigenvalues. 

Although the conjugate pair of eigenvalues lies outside 
the control bandwidth, their inclusion markedly improves the 
output correspondence with respect to the full order system, 


see Table III-2. Further analysis will use the third order 


TABLE III-3 


Reduced Model Eigenvalues 


Model #1 Model #3 Model #2 
XL x1 Xl 
States 
x2 x2 X2 
x5 x5 XS 
x10 X11 
X13 x16 
-7.175 -51.65 + -6.947 
Eigen- 
-3.392 j7.14 -3.282 
Values 
-1.418 - 1.388 -1.389 
- .651 - 6.909 
-1.960 - 3.265 
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Fifth Order Model 
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model (model #2) and the fifth order augmented model (model 


#3). These models are listed in Tables III-4 and III-~S. 


E. VERIFICATION OF REDUCED ORDER MODELS 

The reduced order model must be sufficiently representa- 
tive of the original system for the designed control to pro- 
vide the desired full system response. Weinberg and Adams 
{9], used the F100 non-linear simulation to develop a linear 
17th order system at the static, sea-level point at partial 
power settings. Their reductions to 5th and 3rd order and 
similar reduction to 5th order by DeHoff [10], from Miller 
and Hackney's 16th order model show excellent correlation to 
the non-linear simulation. 

The operating point chosen for consideration in this 
thesis is 30,000 feet, Mach = 0.9 and power level angle of 
67 degrees. This corresponds to a near optimum cruise point 
and is a condition frequently encountered in routine opera- 
tions. Section III.D lists the states and associated eigen- 
values at this operating point and discussed several reduced 
order models. The dominant eigenvalues were considered to 
be those which affected the models fast response and had 
dynamic significance to the engine's physical characteristics. 
Again, this implies a control bandwidth of 0-10 rad/sec, how- 
ever, other models, using eigenvalues outside the control 
bandwidth, were considered in determining the best reduced 
model for deriving the control laws. 

One method to analyze model correspondence is to use the 


Bode plot across the control bandwidth. Figure III-l is a 
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Bode comparison of the DeHoff reduced model and the Miller 
and Hackney 16th order model at static, sea-level, inter- 
‘ ; mediate p47r for the 0 to 10 rad/sec control bandwidth. 
It is clear that through the 0-10 Hertz bandwidth the two 
models show very close correspondence; this similarity breaks 
down in the higher frequency region, but the response time of 
the eigenvalues outside the control bandwidth is considered ‘ 


to be too short to affect transient behavior. The phase 


response, not shown, follows the same trends. 


Iv. ANALYSIS OF CONTROL SYSTEMS 


A. COST FUNCTION WEIGHING FACTORS 


The control law for the engine regulating systen, 


U = Un + C,, (X-X) + | CYA (Y“Yo) at (IV-1) 
can be reduced as shown in Section III.B to 


U = U_ = KX (IV-2) 


when considering small perturbations and ignoring the out- 

puts from the system. In most circumstances, the Un! the 

trim condition, would be considered as zero for steady state 
analysis. The linear models developed by Miller and Hackney 
[5], set all inputs to zero at the steady-state operating point 
for which the model is defined. 

The K matrix is the solution to the optimum regulator 
problem discussed in Section III.A. The OPTSYS4 computer 
solution used in this thesis uses the cost function 
cc; 


ae er : 
Jo= Z/(Y AY + UR Uldt (IV=3) 


where A_ is the weighting matrix for the output vector and 
| replaces the need for the Ry Matrix in equation (III-7). 
Kirk {1ll], discusses the effect that weighting matrices 


have on system performance. In general, the larger the 
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Magnitude of the matrix element the faster its associated 
parameter (system state, input or output) stabilizes. An 

’ additional factor in determining the magnitude of the diagonal 
Matrix elements is the magnitude of their associated param- 


eters. Given the cost function 


= l T. J 
J = x f X°R)X dt (IV=4) 
where 
1 
0 xl 
ae 
Ry = : X = (IV=-5) 
9) + x2 
b 


The cost function then becomes 


X2,2 
soy } dat (IV-6) 


= xl,2 
J = {( _ + ( 
The desire is to minimize the cost function, J, with 
respect to time. Thus, taking the time derivative and 
setting it equal to zero 


+ = 0 (IV=7) 


xl,2 X2, 2 
=) (>) 
This is the solution which produces a minimum cost. The 


values of a and b are the largest expected values of xl and 


ng 


x2, respectively, and represent an equal effect on the cost 
function by each state. Should either a or b be decreased, 
: . the optimal control developed would have to increase the 


control input to the other state to minimize the cost func- 


tion. A nominal fluctuation limit of +5% will be considered 

as a small perturbation in the initial control analysis. 
Designing a control to minimize the thrust specific fuel 

consumption (SFC) of the engine requires limiting the thrust 

and the fuel flow fluctuations. The weighting matrix ele- 

ments corresponding to these two parameters are, therefore, 

to be increased with respect to the other parameters. Looking 

at only the small perturbation response, it is noted that 

very little change in the total engine airflow (y3) is ex- 

pected and, as before, there is no change in customer bleed 

air requirements (u5). These two parameters can be neglected, 

as their response characteristics have no effect on the desired 

SFC response. To demonstrate the effect of the different ele- 

ment magnitudes on the control gain matrix, the OPTSYS4 pro- 

gram was run for the 3 cases delineated in Table IV-l. Case 

I was developed to allow each output and input, with the excep- 

tion of the two pressure ratios (y6 and y7) and the bleed air 


percentage (u5), to effect equal contribution to the cost 


7 function. The diagonal elements of both ay and Ro normalize 
the squared value to the order of one. Case II uses all of 
the same elements, except the thrust (yl) and the fuel flow 


(ul) are weighted more heavily, to the magnitude of ten. This 


TABLE Iv-1 
Case Definition 


A. and R, are diagonal matrices, only the significant 
elements are listed. 


Case I 
AY = (.0001 1 .0001 100 100 le-8 le=8) 
Ro = (.0001 100 .01 -01 le-8) 

Case II 
Ay = (.01 1 .0001 100 100 le=-8 Ile-8) 
Ro = (.01 100 .01 -Ol1 Le~8) 


Case III 
A = (1 1 .0001 100 100 tIle=8 le-8) 
a ?0'- 00. 08 -01 le=8) 


heavier weighting should reduce the fuel input to the system 
and moderate the thrust fluctuations as the control returns 
the system to the steady state condition. Case III increases 
the weighting of yl and ul to achieve a normalized magnitude 
of 100. This extremely heavy weighting should eliminate fuel 
from the control and reduce thrust variations to absolute 
minimum. The whole intention of the controller design is to 
damp out state perturbations without using fuel or altering 
the thrust level, i.e., making no change in specific fuel 


consumption to effect control. 
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These three cases were implemented on the two reduced 


order models and simulations run for a perturbation of -50 
RPM in fan speed (xl). The fifth order model returned the 
fan speed to 95% of steady state in .5 seconds using Case I, 
but increased fuel flow by 10 pounds per hour initially and 
reguired an increase in fuel consumption of .2% over the 

half second period. For Case II, the fifth order model re- 
turned to 95% of the steady state fan speed in approximately 
-3 seconds, while using almost 80% less fuel than the Case I 
controller. This 80% reduction in fuel consumption indicates 
that the increased weighting of Case II has accomplished 
exactly what was desired. Case III was run for the fifth 
order model and, although virtually no fuel was used in the 
control, the system experienced oscillations in fan speed of 

+8 RPM through the 3 second simulation run. Additionally, 

the Case III controller caused rapid changes of +4 degrees in 
the compressor guide vanes (u4), the regulatory limit imposed 
for the simulation, and also large fluctuations in the inlet 
guide vanes (u3). One would achieve results of similar magni- 
tude for each of the cases assuming a +1% perturbation of all 
states in the reduced-order model. Case I is too fuel de- 
pendent and causes large changes in u3, Case III also requires 
large changes in U3 that could cause a slew rate request larger 
than that available from the control actuator. 


In modeling the control, the actuator and physical limi- 


tations of each control must be considered and appropriate 


limitations placed on the control law to avoid overshoots. 


Table IV-2 gives the actuator limits and maximum rates 


applicable to the F100 engine [4]. A stepping limit of .01 


seconds is assumed in determining the rate limits and fuel 


flow fluctuations are limited to +10% of the steady-state 


value. No compensation is allowed for hysteresis or non- 


linear dynamics of the actuators; these are considered 


minimal in the small perturbation case, but would need to 


be incorporated for large perturbation control modeling. 


TABLE IV-2 
Actuator Limitations 


Maximum Minimum Rate Limit 


Input 
Fuel flow 16300 lb/hr 450 lb/hr 15800 lb/hr/sec 


Nozzle area 6.4 sq. ft. 2.8 sq. ft. 3.6 sq.ft./sec 


-40 deg. 48 deg/sec 


Inlet guide vanes 0 deg. 


-40 deg. 40 deg/sec 


4 deg. 


Compressor vanes 


It is evident that the Case III weighting is not accepta- 


large amplitude 


ble because of the oscillations and the rapid, 


fluctuations in the guide vane positions (u3 and u4). Casel 


is also inappropriate for the specific fuel consumption mini- 


mizing control since it uses fuel as the primary control in- 


put. For these reasons, Case II is chosen as the weighting 


for all control modeling to be implemented in later sections. 


and the feed- 


Once the integral portion of equation (IV-1) 


back gains are determined, the constant portion, Un! must be 


analyzed. In most situations this term would be neglected. 
The steady state analysis is now performed to determine if 
the specific fuel consumption can be reduced by making small 
changes to the trim position of the variable geometry. An 
analysis of control changes versus SFC, as was done at the 
static sea-level condition in Section II.D, must again be 


conducted. A minimum SFC obtainable relates to certain input 


values in u. which is incorporated into the control law, im- 


proving the steady-state performance and controlling the 
fluctuations of the system which will have the greatest effect 
on the SFC. 

As was done at the static sea-level operating point, a 
small perturbation in each input was made and its effect on 
the SFC, thrust and other engine conditions tabulated. The 


goal of these analyses is to develop a trim input, U that 


Im? 
improves the fuel economy at steady state without causing 
significant change in the operating conditions (both states 
and outputs). 

The effects of the control inputs on several parameters 
are plotted on Figures IV-l through Iv-4. Figures IV-1l 
shows, again, the large effect nozzle area (u2) and fuel flow 
(ul) have on specific fuel consumption, the effect of inlet 
guide vanes (u3) is drastically reduced from that observed 
at static, sea level operation. Figure IV-2 indicates the 


offsetting effect of ul and u2 and again indicates a zero 


net thrust control can be found. Fan speed (xl) and compressor 
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speed (x2) effects are shown in Figures IV-3 and Iv-4, respec- 
tively. At this operating point the inlet guide vane posi- 
tion (u3) may be most useful in controlling the fan speed 
fluctuation and thereby stabilizing the fan stability margin 
(y4). A sample control would close u2 by .2 square feet and 
reduce fuel flow by 57 pounds per hour. This combination 
produces a 1.4% decrease in specific fuel consumption and 
results in a 4% reduction in fan speed. A trim condition 


control ing t, Un would be 


B. REDUCED ORDER SIMULATION ANALYSIS 

The model developed in Section III.c and Section III.D for 
the high altitude flight operating point was used to initially 
analyze the control law developed in Section IV.A. A reduced 
order controller is simpler and cheaper to implement through 
a digital computer. It requires less computer space and far 
less computing time than the full-order linear model developed 
controller. 

A comparison of the three cases for which gain matrices 
were found indicated that Case 2 gave the most desirable re- 
sults. The Case 2 weighting matrix was applied to the third 


order reduced system and perturbations applied to various 
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states. The same weighting matrix was applied to the fifth 


order model (#3) and similar perturbations applied. Figure 
IvV-5 shows the system response to a 50 RPM perturbation of 
the fan speed. The fifth order system has a settling time 
of approximately .45 seconds whereas the third order system 
has one of 1.5 seconds. This variation is accounted for by 
the two fast response poles in the fifth order model, 
(-51.65,37.14). 

These results must be compared to the 16th order system 
to determine which model most closely approximates the full 


order system response. 


C. FULL ORDER SIMULATION ANALYSIS 

The proof of the viability of a reduced order controller 
is its ability to control the full order system in the same 
Manner as predicted in she reduced order model. The two 


reduced mocels analyzed in Section IV.B were applied to the 


sixteenth order iinear model. Additionally, the Case 2 weight- 


ing matrix was applied and a full order controller was found 
and implemented. All simulation was done using the IBM Con- 
tinuous System Modelling Program, Model III. 

All three models were compared with similar perturbations 
applied. Figure IV-6 compares the full order and fifth order 
response to the 50 RPM speed perturbation. Figure IV-7 com- 
pares the full order and third order response to the same 
perturbation. Figure IV-8 compares the two reduced order 


controllers and the full order controller for the same fan 


speed perturbation. 


asuodsay [TO1qUOD JapPAO TINd SA A®PIO YYFTA °9-AT BINbtd 
(09S) ITLL 


£ 4 T . 
a 5+ 


IApIO UOT 


(WRI) 
paads 
ued7 


ee ae 


69 


asuodsey [O1}U0D APIO TINA SA ASpAO PATUL “L-AI erzNhbT4 


(90S) SUTL 


€ Z T 0 
lo ee ee en et ee en ff nen nh cee meg OG= 


| 


Pa coe me pee ep ee eee eee be en cm tte OG 


70 


ne ee ee 


uoT}eQUsWEeTdUT JepPIO TInd °*8-AI saznbTd 


(99S) SuUTL 
4 
: a es T 
4 4-- oe es Ps es 0S 
eee BE sere Looe 0 janis = 
ue dV 


The results show satisfactory agreement between the re- 
duced order and full order implementation of the reduced 
order controller. This proves the viability of the con- 
troller and implies that comparable results will be found 
if the control is implemented in the non-linear dynamic 
syaten: The fifth order controller, with two poles outside 
the control bandwidth, most closely models the response of 
the full order controller. As was seen previously, the 
inclusion of the conjugate pair of poles dramatically im- 
proves the response time of the model and provides a much 
closer approximation to the full order system Pepeounen 

This indicates that the fifth order model is sufficient 
for the control development process at this operating point. 
The two fast response poles outside the control bandwidth 


must be included to achieve this close simulation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The study of implementing a reduced order controller 
that is flight regime dependent shows that at a specific 
operating point an improvement in specific fuel consumption 
can be made. The conclusions reached at the end of this 
research are: 

1. The variable geometry existing in the state-of-the- 
art turbofan engines can be used to further improve the 
specific fuel consumption; 


2. A reduced order controller can be found and imple- 


mented that concentrates control authority on the variable 
geometry rather than the fuel control, thus furthering the 
fuel savings; and 

3. A similar controller can be found and implemented at 
any given operating point throughout the flight regime, these 
controllers could then be organized into a single digital 
machine to provide a regime-dependent engine control to im- 


' prove the airframe/engine interaction in all flight conditions. 


B. RECOMMENDATIONS FOR FUTURE STUDY 

One obvious extension of this work is to find further 
reduced order controllers. 

Additionally, work can be done to fully implement these 
reduced order controllers in the non-linear model using 


either the Pratt & Whitney transient simulation digital model 


or the hybrid computer model. 


The effect on fatigue life and life cycle cost due to 
the "fuel saving" control schemes presented was not deter- 


mined. It is therefore recommended that these factors he 


explored in order to fully analyze the usefulness of these 


control schemes. 
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